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The vibron dynamics associated to amide-I vibrations in a 3D a-helix is described according to a 
generalized Davydov model. The helix is modeled by three spines of hydrogen-bonded peptide units 
linked via covalent bonds. To remove the intramolecular anharmonicity of each amide-I mode and 
to renormalized the vibron-phonon coupling, two unitary transformation have been applied to reach 
the dressed anharmonic vibron point of view. It is shown that the vibron dynamics results from 
the competition between inter-spine and intra-spine vibron hops and that the two kinds of hopping 
processes do not experience the same dressing mechanism. Therefore, at low temperature (or weak 
vibron-phonon coupling) , the polaron behaves as an undressed vibron delocalized over all the spines 
whereas at biological temperature (or strong vibron-phonon coupling), the dressing effect strongly 
reduces the vibrational exchanges between different spines. As a result the polaron propagates along 
a single spine as in the ID Davydov model. Although the helix supports both acoustical and optical 
phonons, this feature originates in the coupling between the vibron and the acoustical phonons, 
only. Finally, the lattice distortion which accompanies the polaron has been determined and it is 
shown that residues located on the excited spine are subjected to a stronger deformation than the 
other residues. 

PACS numbers: 63.20.Dj,63.22.+m,71.38.Ht,87.10.+c 



I. INTRODUCTION 

In living systems, the energy released by the hydroly- 
sis of adenosine triphophate (ATP) is a universal energy 
source allowing many biological processes such as muscle 
contraction, active transport or enzyme catalysis. How- 
ever, the fundamental question arises whether this energy 
can be transported from the active sites of the living cell 
to other regions without being dispersed or dissipated. 

This question was first pointed out by Davydov and 
co-workers in the 70th to explain the energy transport 
in a-helices 0, 0- The main idea is that the released 
energy, partially stored in the high-frequency amide-I vi- 
bration of a peptide group, can be transported from one 
end of the helix to the other. The dipole-dipole coupling 
between the different peptide groups leads to the der- 
ealization of the internal vibrations and to the formation 
of vibrons. However, the interaction between the vibrons 
and the phonons of the helix induces a nonlinear dynam- 
ics which counterbalances the dispersion and yields the 
creation of the so-called Davydov soliton (for a recent 
review, see for instance Refs. UQ). 

Unfortunately, no clear evidence has yet been found 
for the existence of solitons in real proteins and it has 
been suggested by Brown and co-workers Q and by 
Ivic and co-workers 0, H, E| that the solution is rather 
a small polaron than a soliton. Indeed, this problem 
exhibits two asymptotic solutions depending on the val- 
ues taken by three relevant parameters, i.e. the vibron 
bandwidth, the phonon cutoff frequency and the strength 
of the vibron phonon coupling. When the vibron band- 
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width is greater than the phonon cutoff frequency, the 
adiabatic limit is reached so that the phonons behave in 
a classical way and create a quasistatic potential well re- 
sponsible for the trapping of the vibron according to the 
Davydov soliton. By contrast, when the vibron band- 
width is lesser than the phonon cutoff frequency, the sit- 
uation corresponds to the non-adiabatic limit in which 
the quantum behavior of the phonons plays a crucial role. 
A vibron is thus dressed by a virtual cloud of phonons 
which yields a lattice distortion essentially located on a 
single site and which follows instantaneously the vibron. 
The vibron dressed by the lattice distortion forms the 
small polaron and, as mentioned by the authors, this sit- 
uation corresponds to protein dynamics. Recently, the 
polaron approach has been improved to characterize the 
two-polaron energy spectrum with a special emphasis 
onto the interplay between the dressing effect and the 
intramolecular anharmonicity of each amide-I vibration 
[IcllTlj . Both features favor the formation of two-polaron 
bound states and it has been shown that the helix sup- 
ports two kinds of bound states referring to the trap- 
ping of two polarons onto the same amide-I mode and 
onto two nearest neighbor amide-I vibrations. These re- 
sults were corroborated by a recent experiment devoted 
to the femtosecond infrared pump-probe spectroscopy of 
the N-H mode in a stable a-helix which has revealed the 
two excited state absorption bands connected to the two 
kinds of bound states 0, ^| . Vibrational self-trapping 
in model helix was thus observed for the first time, vali- 
dating in the same time the polaron approach. 

However, most of the theoretical studies applied to the 
Davydov problem involve a ID approximation of the 3D 
structure of the a-helix. Although a real helix is formed 
by three spines of hydrogen-bonded peptide units con- 
nected through covalent bonds, such an assumption was 
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motivated by the fact that when a vibron is symmet- 
rically excited over the three spines, the vibron-phonon 
dynamics reduces to that of a single hydrogen bonded 
chain 0,0, 0,^1- However, it has been shown within 
the soliton approach, that when a vibron is excited on 
a single spine the 3D nature of the helix manifests itself 
by reducing the soliton velocity and by favoring the ex- 
changes between spines [lj, [15J, [l6| . Recently, a detailed 
analysis of both the soliton trapping and the soliton dy- 
namics in a 3D a-helix has been done by Hennig jl7j . 
It has been shown that the soliton envelop exhibits a 
monotonic decay but shows a multihump structure. In 
addition, the lattice distortion along hydrogen bonds is 
reduced for the 3D model when compared with the ID 
model, and, in the excited region, the radius of the helix 
decreases. From a dynamical point of view, Hennig has 
recovered the soliton velocity decrease and he has shown 
that the preferred transport pathway takes place along 
the hydrogen bonds. 

In this paper, the small polaron formalism is applied 
to a 3D model of an a-helix with a special emphasis onto 
the modification of the dressing mechanism induced by 
the 3D nature of the phonons of the helix. We restrict 
our attention to the characterization of the single-polaron 
energy spectrum and the generalization to two-polaron 
eige nstates is addressed in the second paper of this series 

The paper is organized as follows. The model to de- 
scribe the vibron-phonon dynamics in 3D a-helices is in- 
troduced in Sec. II. In Sec. Ill, the nature of the phonons 
in the helix is first summarized. Then, we first realize a 
unitary transformation to remove the intramolecular an- 
harmonicity and a Lang-Firsov |19| transformation is ap- 
plied to renormalize the vibron-phonon interaction. Fi- 
nally, a mean field procedure is performed to obtain the 
effective dressed vibron Hamiltonian. In Sec. IV, a nu- 
merical analysis of the vibron-phonon dynamics is per- 
formed depending on the values of the relevant param- 
eters of the problem. Finally, the results are discussed 
and interpreted in Sec. V. 



z/h 



(a) 
14 

12 

10 

8 
6 
4 
2 




A'' 



A 



A 



A 



A 



A 



A 




-1 -0.5 0.5 1 

x/R 

FIG. 1: (a) Representation of the a-helix backbone denned 
by the angle 0o = 100°, the translation h = 1.5 A and the 
radius Ro = 2.8 A. Each site n is occupied by a residue and 
the dashed line characterizes a particular spine of hydrogen- 
bonded peptide units, (b) Top view of the helix backbone 
where a particular local frame is displayed. 



II. HELIX STRUCTURE AND MODEL 
HAMILTONIAN 

Let us consider a sequence of N amino acid units 
(called residues) regularly distributed along a polypep- 
tide chain. Let n the index which labels the position 
of the nth residue along the chain. In three dimension, 
the structure of the polypeptide chain is stabilized by 
the hydrogen bonds between the carboxyl oxygen (CO) 
of residue n and the amide hydrogen (NH) of residue 
n + 3. The resulting conformation is a 3D a-helix in 
which each residue is related to the next one by a trans- 
lation of ft. = 1.5 A and by a rotation of 9q — 100° leading 
to 3.6 residues per turn of helix (see Fig. la). The pitch 
of the helix, i.e. the product of the translation (1.5 A) 
by the number of residues per turn (3.6), is equal to 5.4 



A. The radius of the helix, which corresponds to the dis- 
tribution of the center of mass of the residues, is fixed to 
Ro = 2.8 A 0. 

To model the vibron-phonon dynamics in a rather sim- 
ple way, we follow the procedure introduced by Hennig 
|17| and treat the N residues as point-like entities which 
the equilibrium positions are located on sites distributed 
along the helix. Therefore, as shown in Fig. la, the 
geometry of the helix is specified in the (xyz) cartesian 
frame where the z direction is parallel to the axis of the 
helix. The equilibrium position of the nth residue is thus 
expressed as 

R(n) = Rq cos n8oe x + Rq sin u8q e y + nhe z ( 1 ) 
where e a denotes the unit vector parallel to the direction 
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a = x, y, z. 

The nth site contains a amide-I vibration which be- 
haves as an internal high frequency oscillator described 
by the standard creation and annihilation vibron opera- 
tors 6+ and b n . By introducing the intramolecular anhar- 
monicity of each amide-I mode, the vibron Hamiltonian 
is written as 

H v = Y, h ^obl l b n + h- f3 (bl l + b n ) 3 + h-f 4 (bl + b n f + ... 

71 

- s EE -»')&+ &»)(&+&»o (2) 

n n' 

where luq stands for the internal frequency of the nth 
amide-I mode and where J(n — n') denotes the lateral 
hopping constant between the residues n and n'. In 
Eq-(0), 73 and 74 represent the cubic and quartic an- 
harmonic parameters of each amide-I mode. Note that 
for the single-vibron dynamics, the intramolecular an- 
harmonicity does not play a crucial role. As it will be 
shown in the following of the text, it is responsible for 
a renormalization of the harmonic parameters (internal 
frequency, hopping constants, vibron-phonon coupling) 
but it does not fundamentally change the physical be- 
havior of the single-vibron states. Nevertheless, it has 
been introduced here because it drastically affects the 
two- vibron dynamics, as shown in paper II [18j. 

The amide-I vibrations interact with the phonons of 
the helix associated to the collective dynamics of the ex- 
ternal motions of the peptide groups. As mentioned by 
Hennig ^tJi both the covalent and the hydrogen bonds 
are assumed to be represented by point-point interac- 
tions. The coupling between two residues is thus de- 
scribed by a pair potential F[r(nn')] which depends on 
the instantaneous distance r(nn') =| r(n) — r(n') | be- 
tween the two residues n and n'. Within the harmonic 
approximation, each peptide group performs small dis- 
placements u(n) around its equilibrium position R(n) so 
that the instantaneous position r(n) of the nth residue 
can be expressed as r(n) = R(n) + u{n). Therefore, by 
expanding the intermolecular potential around the equi- 
librium helix conformation, the phonon Hamiltonian is 
written as 

H p = Y,^ + 1 E $<*fs(nn')u a (n)u (n') (3) 

na na,n' 

where M denotes the average mass of each residue and 
where p(n) is the momentum associated to the displace- 
ment it(n). Both a and /3, equal to x, y or z, account for 
the representation of the vectors in the cartesian frame 
(x,y,z). In Eq.Q, #(nn') is the force constant tensor 
defined as the second derivative of the potential at equi- 
librium and written as 

_ , d 2 V[r(nn')} 
ou a (n)oup{n'j 

According to the Davydov model, the vibron-phonon 
interaction originates in the modulation of the vibra- 



tional frequency of each amide-I vibration by the exter- 
nal motion of the residues. This coupling is expressed in 
terms of the instantaneous positions r(n) of the residues 
as 

AH »p = E \ S K{r{nn')}){bl + b n f (5) 

n 

Within the deformation potential approximation, the 
function Sh({r(nn')}) can be expanded around the equi- 
librium position of the residues so that the resulting cou- 
pling Hamiltonian is finally expressed as 

Ag, = E Xl T l R(nn ltln Un ' ] (6 - + bnf (6) 

n,n' 

where R(nn') — R(n) — R(n') and where X\n—n'\ denotes 
the strength of the interaction between the nth amide-I 
mode and the external displacement of the n'th residue. 

Finally, the vibron-phonon dynamics is described by 
the full Hamiltonian H = H v +H p + AH vp which appears 
as a generalization of the well-known Davydov Hamilto- 
nian to a 3D a-helix. 



III. THEORETICAL BACKGROUND 

In helices, the vibron-phonon interaction is assumed 
to be strong enough so that the vibron dynamics is es- 
sentially governed by the so-called dressing effect. The 
characterization of this effect requires the knowledge of 
the phonon eigenstates and of the Lang-Firsov transfor- 
mation ^3 which yields the small polaron point of view. 
However, although these two features are well known in 
an infinite lattice with translational invariance, the heli- 
cal nature of the lattice is responsible for a modification 
of both the phonon structure and the dressing mecha- 
nism. The present section is thus devoted to the presen- 
tation of these modifications. 



A. Phonons in a 3D a-helix 

In a recent paper devoted to solitons in an isolated 
helix chain, Christiansen et al. have introduced a gen- 
eral formalism to characterize the phonon dynamics in 
an a-helix J2jj ■ This formalism consists in expressing the 
residue dynamics in a simple way by taking advantage of 
the invariance of the phonon Hamiltonian according to a 
rotation by an angle of 9q around the helix axis followed 
by a translation along this axis by a distance h plf. 

To proceed, it is more convenient to express the 
position of the residues in a set of N local frames 
(e r (n), eg(n), e z (n)) attached to each residue (see Fig. 
lb). The unit vectors of the nth local frame are ex- 
pressed in terms of the unit vectors of the cartesian frame 

{&xi t~"yi ^2) 

e M (n) = l) IQ (n)e Q (7) 
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where fi — r,8, z and a — x, y, z and where T(n) is the 
local (3 x 3) rotational matrix defined as 



cosn#o sin n(9 0" 1 
T(n) = | — sinn# cosn8 
1, 



(8) 



In that context, let v{n) denotes the displacement of 
the nth residue expressed in the nth local frame. This 
vector is defined in terms of its cartesian counterpart 
u(n) as v(n) = T(n)u(n). Therefore, by applying this 
latter transformation, the phonon Hamiltonian Eq.Q is 
rewritten as 



H,, 



2M 



(9) 



where p^(n) denotes the momentum connected to v^(n) 
and where *(rm') = T(n)^(nn')T^ 1 (n'). 

At this step, it is straightforward to show that the 
transformed force constant tensor \&(nn') depends only 
on the index difference n — n' (see Appendix A). As a 
consequence, the Bloch theorem can be applied so that 
the displacement of the nth residue in the nth local frame 
can be expanded as a superimposition of plane waves as 



(10) 



where k denotes a reduced wave vector which belongs to 
the first Brillouin zone of the helix, i.e. — ir < k < tt . The 
phonon Hamiltonian is thus defined in terms of the (3 x 3) 
dynamical matrix D^ u {k) = J2 n ^^i® 71 ) exp(ikn)/M as 



k (J, 



2M 



M x ^ 

kfi : v 



D^(k)v^(k)v v (k) (11) 



where p(k) is the momentum associated to v(k). Note 
that the expression of the Dynamical matrix is detailed 
in Appendix A. 

Finally, the normal mode decomposition is achieved by 
performing the diagonalization of the dynamical matrix 
D(k) for each k value. Such a procedure allows us to 
define three eigenvalues f^ CT and three eigenvectors e ka 
labeled by the index a — 1,2,3. Therefore, the two in- 
dexes k and a specify a particular phonon mode with 
energy M7fc CT , quasi- momentum hk and polarization e ka . 
The quantum dynamics of each mode is described by the 
well-known creation a\ a and annihilation a k(J operators 
so that the phonon Hamiltonian can be rewritten in the 
standard form as 



ka 



(12) 



In the same way, the displacement of the nth residue is 
finally expressed as 



*v( n ) 



£ 

ka 



2MNVt ka 



(a ko 



-kcr ' 



^■fika^ 



(13) 



At this step, the vibron-phonon coupling Hamiltonian 
Eq.© can be rewritten in terms of the phonon normal 
mode coordinates as 



A 



ka 



^e* kn a ka ){bi + b n f 



where Afc CT accounts for the modulation of the frequency 
of the nth amide-I vibration due to its coupling with 
the phonon mode specified by the wave vector k and the 
index a. It is defined as 



X\ m \C k {m)e* ko 



^J2MNhtt k<J 
where Cfc(m) denotes a vector expressed as 



Cfc(n'-n) = [T(n) - T(n')e* 



R{nri) 



(15) 



(16) 



B. 



Anharmonic vibron, small polaron and effective 
Hamiltonian 



To remove the cubic and quartic intramolecular anhar- 
monicity, the standard procedure described by Kimball 
et al. |22| is used. First, by disregarding the lateral cou- 
pling between nearest neighbor amide-I vibrations as well 
as the vibron-phonon interaction, the procedure consists 
in performing a unitary transformation V to diagonalize 
each anharmonic amide-I mode. Then, an approximate 
Hamiltonian is obtained by applying the transformation 
on the full Hamiltonian H and be keeping the vibron- 
conserving terms, only. As a result, by using the proce- 
dure detailed in Ref. 0, , the transformed Hamilto- 
nian H = VHV + is expressed as (within the convention 

h=i) 

H = J2^o-2A-B)bl l b n ~J2 J i( n - n ') b n b n'+H P 

n nn' 

+ ^(l + 2n)(A to e- tt "aL + ALe lfc "a feCT )6t6„ (17) 
No- 
where A = 3O7I/W0 — 674 denotes the positive anhar- 
monic parameter and where the other parameters in 
Eci. (|17|) are defined as 



B = J2 72 J (n) (73/ 



u ) 



Jx(n) = J(n)(l + 44(7 3 /^ ) 2 - 12 7 4/^ ) 

n = 120(73/^)2 ~ 12 74 M, (18) 

Note that since we restrict our attention to the single- 
vibron dynamics, additional terms which involve the 
product of four vibron operators have been disregarded 
in the present version of the Hamiltonian Eq. (|17|) . These 
terms play a crucial role for the two-vibron dynamics 
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[13 E3 and they are addressed in the second paper of 
this series [l8j . 

To partially remove the vibron-phonon coupling 
Hamiltonian, a Lang-Firsov transformation is applied 
|19|. Indeed, since the vibron-phonon dynamics is domi- 
nated by the so-called dressing effect, we consider a "full 
dressing" and introduce the following unitary transfor- 
mation 



(7 = exp(^(l + 2ry)[ 



Afc CT e 



-ikn 



A? ft 



ikn 



Qka ka Qka 



-a-k 



(19) 



By using Eq. l|19() . the transformed Hamiltonian H = 
U HW is written as 

H = H p + J2^o-2A-B-(l + 4r])E B )bib n 

Th 

■ Y,Mn-n')&le n ,btb n , (20) 



n') exp(— (1 + Ari)S(n — n',T)) is expressed in terms of 
the coupling constant S(n,T) defined as (ks denotes the 
Boltzmann constant) 



S(n,T) = X] 



ka 



A ka ,2 
Qka 



coth(^%)(l - cos(fcn)) (25) 
2k B ± 



The Hamiltonian -ff e // (Eq.lPU) describes the dynam- 
ics of anharmonic vibrons dressed by a virtual cloud of 
phonons, i.e. small polarons. A polaron state corre- 
sponds to a delocalized plane wave with wave vector q 
and frequency tu q — a> — 2 J2 n >i Jeff( n ) cos(qn). Eq.J53J 
accounts for a renormalization of the main part of the 
vibron-phonon coupling within the non-adiabatic limit. 
The remaining vibron-phonon coupling is thus assumed 
to be small in order to be treated using perturbative the- 
ory. Nevertheless, such a contribution is disregarded in 
the present work. 



where E B denotes the small polaron binding energy de- 
fined as 



(21) 



whereas 8„ stands for the dressing operator expressed as 



e, 



exp(- 



(l + 2r,)[ 



A fc , 



feo-e 



-ikn 



A* P ikn 



n 



ka 



l ka 



ka 



Oka]) 

(22) 

In this dressed vibron point of view (Ea.l|2()|l). the 
vibron-phonon coupling remains through the modulation 
of the lateral terms by the dressing operators. Although 
these operators depend on the phonon coordinates in a 
highly nonlinear way, the vibron-phonon interaction has 
been strongly reduced within this transformation. As a 
result, we can take advantage of such a reduction to per- 
form a mean field procedure 7, 8, 23] and to express the 
full Hamiltonian H as the sum of three separated contri- 
butions as 



H = H P 



IS 



H n + AH 



(23) 



where H e ff = ((H — H p )) denotes the effective Hamilto- 
nian of the dressed vibrons and where AH = H — H p — 
((H — H p )) stands for the remaining part of the vibron- 
phonon interaction. The symbol (. . .) represents a ther- 
mal average over the phonon degrees of freedom which 
are assumed to be in thermal equilibrium at temperature 
T. 

As a result, the effective dressed vibron Hamiltonian 
is written as 



where (2>o — uiq — 2 A — B — (1 + 4t])Eb and where the 
effective vibron hopping constant J e f f(n — n') = J\(n — 



IV. NUMERICAL RESULTS 

In this section, the previous formalism is applied to 
characterize the single-vibron energy spectrum of a 3D 
a-helix. However, the present theory involves a set of 
parameters which has first to be discussed. 

From the literature, the parameters enter the vibron 
dynamics in 3D a-helices are relatively well described . 
The quantum energy for an amide-I vibration is about 
1660 cm -1 so that the harmonic frequency is fixed to 
ujq = 1695 cm -1 . A well admitted value for the hop- 
ping constant along the hydrogen bonds is J(3) = 7.8 



cm (sec for instance Rcfs. 



]). In the same way, 



the vibron hopping constants between different spines of 
hydrogen-bonded peptide units have been calculated in 
Ref. 24] so that J(l) = -12.4 cm" 1 and J(2) = 3.9 
cm -1 . Despite their small values, the other elements of 
the hopping matrix have been taken into account in the 
present work and are listed in the Table I of Ref. [24| . Fi- 
nally, the anharmonic constant of the amide-I vibration 
is about A = 8.0 cm" 1 [13, |M and both the cu- 
bic and quartic anharmonic parameters can be expressed 
approximately by using the relation 15j^/ujq « 674 ~ A 
[T(j. Note that these parameter values are in a rather 
good agreement with recent ab initio calculations |27J. 

The phonon dynamics is essentially governed by the 
force constant tensor defined in Eq.QJ. In that context, 
we assume that the helical conformation corresponds to 
the equilibrium of each pair potential which connects 
two residues. Therefore, the force constant tensor is ex- 
pressed as 



= -K\r 



R a (nn')Rp(nn') 
R 2 (nn') 



(26) 



where K\ n _ n i\ denotes the second derivative of the pair 
potential connecting the two residues located onto the 
sites n and n'. According to the 3D model of Hennig 
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FIG. 2: Phonon dispersion curves for Ki = 60 Nm" 1 , K 3 = 
15 Nm" 1 (full line) and for Ki = 75 Nm" 1 , K 3 = 20 Nm" 1 
(dashed line) (Note that K2 = 0). The dotted line represents 
the phonon dispersion curve of the ID Davydov model with 
Ks = 15 Nm -1 . 



|17|. the a-helix exhibits two modes corresponding to a 
longitudinal displacement along the helix axis and to a 
radial displacement perpendicular to this axis. To re- 
cover these results, the force constant A2 is set to zero 
and we assume K\ n \ =0 for n > 4. The force constant A3 
of the hydrogen bonds ranges between 13 and 20 Nm -1 
Pi IIU fill fl7| whereas the force constant K\ of the co- 
valent bonds varies between 45 and 75 Nm -1 [13. With 
our notations, K\ as W/2, where W is the force constant 
for the radial motion introduced by Hennig. Finally, the 
mass M which enters the phonon dynamics has been fixed 
to 2.0 10~ 25 kg. 

The vibron-phonon coupling is basically defined from 
the knowledge of the parameters X\n-n'\ (Eq-©) which 
accounts for the modulation of the nth amide-I frequency 
due to the external motion of the n'th residue. Accord- 
ing to the ID Davydov model, the parameter ^3 ranges 
between 35 and 62 pN. However, the influence of the mo- 
tion of the residues n ± 1 and n ± 2 onto the nth amide-I 
vibration is still unknown. Therefore, we choose X2 = 
whereas we treat %i as a parameter smaller than ^3 |l7j . 
Note that \n = for n > 4. 

In Fig. 2, the phonon energy spectrum is shown for 
A'i = 60 Nm -1 , A 3 = 15 Nm" 1 (full line) and for 
A'i = 75 Nm -1 , A 3 = 20 Nm" 1 (dashed line). The 
spectrum exhibits two non vanishing dispersion curves 
whereas the third one is equal to zero in the entire Bril- 
louin zone. The low frequency curve, which tends to zero 
when k — ► 0, corresponds to acoustical phonons, whereas 
the high frequency curve represents optical phonons. In 
the long- wavelength limit, the acoustic branch increases 
linearly with the wave vector k and describes a sound 
wave. However, the curve reaches a maximum value when 



k = 7r/3 then it decreases to a minimum value when 
k = 2tt/3. Finally, the frequency increases again when 
k reaches the edge of the Brillouin zone. Note that the 
acoustic branch is rather insensitive to the force constant 
Ai whereas it strongly depends on the force constant A3. 
The maximum acoustic frequencies are equal to 89 and 
102 cm -1 when A3 = 15 and 20 Nm -1 , respectively, 
whereas the minimum frequency in k = 2tt/3 is about 
20 cm -1 . The optical branch appears almost constant 
over the entire Brillouin zone and exhibits a rather weak 
dispersion of about 8.5 cm -1 . It strongly depends on 
the force constant Ai and appears almost insensitive to 
A3. At the center of the Brillouin zone, the optical fre- 
quency is equal to 135 cm -1 when Ai = 60 Nm -1 and 
reaches 151 cm -1 when Ai = 75 Nm -1 . Note that the 
dotted line represents the phonon dispersion curve in the 
corresponding ID Davydov model for A3 = 15 Nm . 

In the long- wavelength limit, the analysis of the eigen- 
vectors of the dynamical matrix reveals that the acoustic 
phonons are essentially polarized along the axis of the 
helix. More precisely, they correspond to longitudinal 
vibrations along the direction specified by the hydrogen 
bonds. By contrast, the optical branch refers to a ra- 
dial motion of the residues which describes a breathing 
mode of the radius of the helix. These two features cor- 
roborate the previous observed dependence of the two 
branches with respect to the force constants Ai and A3. 
Note that a singular behavior takes place when k — 2ir/3 
since the acoustic branch becomes polarized along the ra- 
dial direction, only. In the same way, when k reaches the 
edge of the Brillouin zone, the optical branch refers to an 
hybridization between a radial and a torsional motion. 

To characterize the vibron-phonon interaction, the k 
dependence of the normalized function y/N \ Aka | de- 
fined in Ea. H15fl is displayed in Figs. 3. The phonon 
force constants are fixed to Ai = 60 Nm -1 and A3 = 15 
Nm -1 and two situations are considered for the strength 
of vibron-phonon coupling, i.e xi — 25 pN, X3 — 50 pN 
(foil line) and xi = 10 pN, X3 — 60 pN (dashed line). As 
shown in Fig. 3a, the coupling with acoustical phonons 
strongly depends on the phonon wave vector. It vanishes 
for three k values defined as fc^ 1 -* = 0, = tt/3 and 
= 7t. The coupling exhibits three maxima located 
at fc'W = 0.41, k'W = 1.71 an d fc'(3) = 2.48. When 
Xi = 25 pN and X3 = 50 pN, the maxima lye around 
23 cm -1 whereas they reach 28 cm -1 when xi = 10 pN 
and %3 = 60 pN. In fact, the coupling with acoustical 
phonon depends significantly on X3 and appears almost 
insensitive to xi- Note that the coupling shows a local 
minimum around fc = 27r/3 which correspond to a zero 
coupling in the ID Davydov model (dotted line). 

As shown in Fig. 3b, the coupling with optical phonons 
basically decreases when the phonon wave vector in- 
creases. It is maximum at the center of the Brillouin 
zone and vanishes when fc = tt. When xi = 25 pN and 
X3 = 50 pN, the maximum is equal to 20 cm -1 whereas 
it reduces to 15 cm -1 when xi — 10 pN and X3 — 60 pN. 
In a marked contrast with the previous situation, the 
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FIG. 3: Normalized function VN | A fctT (n) | vs k for K ± = 60 
Nm" 1 , K 3 = 15 Nm" 1 and for X i = 25 pN, X 3 = 50 pN (full 
line) and Xi = 10 pN, X 3 — 60 pN (dashed line), (a) coupling 
with acoustical phonons, (b) coupling with optical phonons. 
The dotted line in Fig. 3a refers to the ID Davydov model. 



coupling with optical phonon depends essentially on X \. 
Note that this coupling shows a local minimum around 
k = tt/3. 

In Figs. 4, the coupling constant S(n,T) defined in 
Eq.lHJ is shown for T = 310 K, (full line), T — 150 K 
(dashed line) and T = 5 K (dotted line) and for K\ = 60 
Nm" 1 , K 3 = 15 Nm- 1 , Xl = 25 pN and X3 = 50 pN. 
The sum over the index a in Eq. I|25|) allows us to express 
the coupling constant as S(n,T) — S a (n,T) + S (n,T) 
where S a (n,T) and S (n,T) denote the contributions of 
the coupling with acoustical and optical phonons, respec- 
tively. 

As illustrated in Fig. 4a, at low temperature, the cou- 
pling constant S a (n,T) appears almost site independent 
and slightly fluctuates around 0.25. When the temper- 
ature increases, the coupling constant increases and ex- 
hibits two distinct behaviors depending on the nature of 
the bond n which experiences the dressing effect. Indeed, 
when n = 3,6,9,..., S a (n,T) accounts for the dressing 
mechanism responsible for a decreases of the vibron hop- 
ping constant between amide-I vibrations along the same 
spine of hydrogen-bonded peptide units. The coupling 
S a (n,T) increases almost linearly with n and leads to a 
screening of the corresponding vibron hopping constants. 



FIG. 4: Coupling constant S(n, T) vs the bond index n for 
T = 310 K, (full line), T — 150 K (dashed line) and T = 5 
K (dotted line). The phonon force constants are fixed to 
K\ = 60 Nm -1 and Kt = 15 Nm -1 and the vibron-phonon 
coupling constants are equal to X i = 25 pN and X3 = 50 pN. 
(a) Contribution of the acoustical phonons, (b) contribution 
of the optical phonons. 



When n = 1,2,4,5,..., S a (n,T) characterizes the dress- 
ing of the hopping constant between amide-I vibrations 
located onto different spines of hydrogen-bonded peptide 
units. As previously, it increases as n increases. The 
figure clearly shows that vibron hops between different 
spines experience a stronger dressing mechanism than the 
hops which take place along the same spine. The differ- 
ence between the two kinds of dressing increases as the 
temperature increases. However, as n increases, the two 
kinds of dressing converge to the same n dependence. In 
other words, long range vibron hops experience the same 
dressing effect whatever the location of the amide-I vi- 
brations involved in the process. 

As shown in Fig. 4b, the contribution of the opti- 
cal phonons S (n,T) is about two orders of magnitude 
smaller than the contribution of acoustical phonons. It 
exhibits a fully different behavior with respect to the 
bond index n since it drastically varies with n until n < 4 
and becomes almost n independent when n > 5. Note 
that as previously, the coupling constant S (n,T) in- 
creases with the temperature. 

The temperature dependence of the effective vibron 
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FIG. 5: Behavior of the effective vibron hopping constants 
| J e ff(n) | for n = 1 (dotted line), n — 2 (dashed line) and 
n = 3 (full line). The parameters used for the calculations 
are K\ — 60 Nm -1 and K3 = 15 Nm -1 . (a) Temperature 
dependence for X \ — 25 pN and X z = 50 pN. (b) Eb depen- 
dence for T — 150 K. In that case, Eb varies according to the 
restriction xi = X3/2- (c) do dependence for T = 150 K. 



hopping constants | J e ff(n) \ for n — 1 (dotted line), 
n = 2 (dashed line) and n = 3 (full line) is drawn in Fig. 
5a. The parameters used for the calculations are K i = 60 
Nm" 1 , K 3 = 15 Nm" 1 , X i = 25 pN and X3 = 50 pN. 
Due to the dressing effect, the three hopping constants 
decrease as the temperature increases. In addition, what- 
ever the temperature, the hopping constant | J e ff{2) | 
is always smaller than the two other hopping constants. 
However, the ratio between | J e //(1) | and | J e //(3) | 
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FIG. 6: Small polaron dispersion curve for T = 310K (full 
line), T = 150.fr (dashed line) and T — 5K (dotted line) and 
for K x = 60 Nm" 1 , K 3 = 15 Nm" 1 , X i = 25 pN and X 3 = 50 
pN. 



strongly depends on the temperature and a transition 
occurs when T = T c = 51 K. When T < T c , the constant 
| J e //(1) | is the dominant contribution. For instance, 
at T = K, J e //(1) = -9.7 cm" 1 , J e //(2) = 3.0 cm" 1 
and J e //(3) = 7.2 cm -1 . By contrast, when T > T c , 
I ^e//(3) I becomes the main hopping constant, and, 
for instance at T = 300 K, J e //(1) = -0.43 cm" 1 , 
J e //(2) = 0.13 cm -1 and J e //(3) = 3.7 cm -1 . A simi- 
lar behavior is illustrated in Fig. 5b for the dependence 
of the hopping constants with respect to the small po- 
laron binding energy Eb- The temperature is fixed to 
T = 150 K and Eb varies according to the restriction 
Xi = Xs/2- Due to the dressing effect, the effective hop- 
ping constants decrease as Eb increases. As previously, 
| J e ff (2) | is smaller than the two other constants and 
a transition occurs when Eb = 3.1 cm -1 which discrim- 
inates between | J e //(1) | and | J e //(3) |. Therefore, 
| J e //(1) | is greater than | J e //(3) | when Eb < 3.1 
cm -1 whereas the opposite feature occurs when Eb > 3.1 
cm -1 . Note that when Eb > 25 cm -1 , both | J e //(1) | 
and | Jeff (2) | almost vanish so that only the hopping 
constant | J e //(3) | remains. Finally, the influence of the 
angle 9o on the vibron hopping constants is illustrated in 
Fig. 5c for T = 150 K. The figure clearly shows that a 
change in 6*o is accompanied by drastic modifications in 
the hopping constants | J e //(1) | and | J e ff(2) \ whereas 
| J e ff (3) | is just slightly modified. When f3 is greater 
than 110°, both J e //(1) and J e //(2) vanish so that only 
vibron hops along a given spine of hydrogen-bond peptide 
units remains. 

The small polaron dispersion curve is displayed in Fig. 
6 for T = 310K (full line), T = 150K (dashed line) and 
T = 5K (dotted line) and for K x = 60 Nm" 1 , K 3 = 15 
Nm -1 , xi = 25 pN and %3 = 50 pN. At the center of 
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the Brillouin zone, the polaron frequency is almost tem- 
perature independent and is located around 1660 cm . 
The dispersion curves exhibit a maximum around q = 1.0 
equal to 1697 cm -1 , 1683 cm" 1 and 1677 cm -1 for T = 5 
K, T = 150 K and T — 310 K, respectively. In the same 
way, they show a minimum around q = 2.15 equal to 
1647 cm" 1 , 1657 cm" 1 and 1661 cm" 1 for T = 5 K, 
T = 150 K and T = 310 K, respectively. The polaron 
bandwidth, defined as the difference between the max- 
imum frequency and the minimum frequency, decreases 
as the temperature increases. It is equal 16 cm" 1 for 
T = 310 and reaches 50 cm" 1 for T = 5 K. Note that the 
polaron frequency exhibits a strong temperature depen- 
dence at the edge of the Brillouin zone. 



As defined in Appendix B (see also Refs. (2J, 29] ), 
the distortion of the helix backbone due to the creation 
of a vibron onto the site no — is displayed in Figs. 7 
for Ki = 60 Nm" 1 , K 3 = 15 Nm" 1 , X i = 25 pN and 
X3 = 50 pN. Fig. 7a clearly shows that the vibron is 
responsible for a decrease of the radius of the helix onto 
the sites no, no ± 3, no ± 6 ... which are linked to each 
other via hydrogen bonds. By contrast, the radius of the 
helix on the sites no±l, no±2, no±4, no±5 ... , slightly 
increases but with smaller amplitudes. Note that the de- 
formation of the radius is exponentially localized around 
the excited site no. The angular distortion of the helix 
which corresponds to a torsional motion of the residues is 
shown in Fig. 7b. Although the site no does not exhibit 
any torsion, oscillations in the deformation of the helix 
backbone take place due to the different kinds of bonds 
between the site no and its neighboring sites. Neverthe- 
less, after the distortion has been propagated over few 
residues, the figure clearly shows that the residues lo- 
cated on the right side of the excited site (n > n ) are 
characterized by an increase of the angle around its 
equilibrium value. By contrast, the opposite feature ap- 
pears in the left side of the excited site (n < no) where 
the 6 angle of each residue is basically reduced from its 
equilibrium value. As a result, a kink-like soliton takes 
place which discriminates between a negative torsion of 
about -0.12 degree and a positive torsion of about +0.12 
degree. Similar features are observed for the deformation 
along the axis of the helix as shown in Fig. 7c. Indeed, 
although the central site does not exhibit any distortion, 
a kink-like soliton takes place. In the left side of the ex- 
cited side, the new equilibrium positions of the residues 
along the helix axis correspond to an increase of the z 
coordinates whereas the opposite feature occurs in the 
right side of the excited site. The curve exhibits oscil- 
lations which clearly show that the distortion essentially 
affects the residues no ± 3, no ± 6, ... linked via hydro- 
gen bonds. Note that the amplitude of the distortion is 
about 0.02 A and appears almost 3.5 times smaller than 
the deformation which occurs in the corresponding ID 
Davydov model (dotted line). 
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FIG. 7: Deformation of the helix backbone due to the cre- 
ation of a vibron onto the site no = for K\ — 60 Nm -1 , 
K 3 = 15 Nm" 1 , xi = 25 pN and X3 = 50 pN. (a) radial de- 
formation, (b) angular torsion and (c) longitudinal distortion. 
The Dotted line in Fig. 7c represents the distortion in the ID 
Davydov model. 



V. DISCUSSION 

In the previous section, the numerical calculations have 
clearly shown that the 3D nature of the a-helix strongly 
modifies the small polaron eigenstates when compared 
with the ID Davydov model. However, the origin of these 
modifications is twofold. First, the vibron dynamics is 
different because an helix allows for vibron hops between 
different spines of hydrogen-bonded peptide units. Then, 
the helix supports both acoustical and optical phonons 
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which participate in a different way to the dressing effect. 

Nevertheless, to discuss and interpret these observed 
features, let us first mention some general remarks about 
the modification of the Davydov problem. Indeed, ac- 
cording to the ID model, the vibron-phonon dynamics 
along a single spine of hydrogen-bonded peptide units 
reduces to that of a chain formed by the sites n, n + 3, 
n + 6 ... . In a similar way, since a real 3D helix is invari- 
ant under a rotation of angle 0q followed by a translation 
of length h, the vibron-phonon dynamics reduces to that 
of a ID lattice formed by all the sites n, n+1, n + 2, n + 3, 
... . As the result, in a 3D helix, both the vibron and 
the phonon reduced wave vectors lye in a Brillouin zone 
three times greater than the Brillouin zone introduced in 
the ID model. In addition, the vibron-phonon dynamics 
is now governed by long range interactions. Interactions 
between nearest and second nearest neighbor sites char- 
acterize coupling between different spines of hydrogen- 
bonded peptide units. By contrast, interactions between 
third nearest neighbor sites account for the coupling in 
a given spine. These latter couplings correspond to the 
nearest neighbor interactions in the ID Davydov model. 

The analysis of the collective motion of the residues has 
revealed that the phonon spectrum exhibits an acoustical 
branch. In the long-wavelength limit, this branch char- 
acterizes the longitudinal sound wave which propagates 
along the spines of hydrogen-bonded peptide units, as 
in the ID Davydov model. However, since the residues 
are allowed to move in a 3D space, new features occur 
in a real a- helix. First, the spectrum exhibits an optical 
branch associated to a high frequency breathing motion 
of the radius of the helix. More precisely, this radial 
motion takes place in the long-wavelength limit whereas 
hybridizations with a torsional motion of the helix back- 
bone occur as when increasing the phonon wave vector. 
In other words, optical phonons originate in the cova- 
lent bonds between different spines of hydrogen-bonded 
peptide units. Then, in a marked contrast with the ID 
Davydov model, the polarization of the acoustic branch 
strongly depends on the phonon wave vector. More pre- 
cisely, a strong hybridization between longitudinal and 
radial motion takes place when k reaches 2ir/3. In that 
case, the acoustical branch does not refer to longitudinal 
sound wave anymore. 

In fact, the wave vector k = 2n/3 plays a key role 
for acoustical phonons. Near k = 2ir/3, the shape of 
the acoustic branch strongly depends on the geometry 
of the helix backbone through the value taken by the 9q 
angle. Indeed, when O = 27r/3, our calculations have 
shown that this branch tends to zero (not drawn in Fig. 
2). This surprising feature can be understood as follows. 
When #o = k = 2n/3, it is straightforward to show that 
the local frames n,n±3,n±6, ... are equivalent and that 
all the motions of the corresponding residues are in phase 
since exp(ifcn) = cxp(ik(n ± 3)) = cxp(ik(n ± 6)) = .... 
Therefore, these motions correspond to a zero frequency 
uniform translation along the axis of the helix so that 
the physics is equivalent to that of the ID model at the 



center of its own Brillouin zone. When 0q decreases to 
reach the value characterizing the a-hclix conformation, 
i.e. 9q = 100°, the equivalence between the local frames 
n, n±3, n±6, ... is broken. As a result, in phase motions 
of the associated residues do not correspond to a uniform 
translation anymore but are connected to an hybridiza- 
tion between radial and torsional motions. Therefore, a 
non vanishing frequency occurs in the spectrum. 

As mentioned previously, the vibron dynamics in a 3D 
a-helix is characterized by two kinds of hopping pro- 
cesses. The first kind refers to vibrational transition be- 
tween different spines of hydrogen-bonded peptide units 
whereas the second kind involves vibron hops along the 
same spine. In that context, the nature of the vibron dy- 
namics essentially results from the competition between 
these two kinds of processes. Within the undressed limit, 
i.e. when the vibron-phonon interaction is disregarded, 
vibron transitions between nearest neighbor and third 
nearest neighbor sites represent the dominant hopping 
mechanism since J(l) = —12.4 cm -1 and J(3) = 7.8 
cm -1 whereas J(2) = 3.9 cm -1 . In other words, the 
3D nature of the helix plays a crucial role and the vi- 
bron is delocalized between the different spines. How- 
ever, our calculations have clearly shown that the dress- 
ing effect strongly modifies the value of the effective hop- 
ping constants and affects in a different way the two 
kinds of hopping processes. Indeed, at low temperature 
| J e //(1) |>| J e ff (3) |>| J e ff (2) | so that the small po- 
laron behaves almost like the undressed vibron. This is 
no longer the case as when the temperature is increased 
because the two kinds of hopping processes do not ex- 
perience the same dressing effect. Indeed, the dressing 
for vibrons hops between different spines is more efficient 
than the dressing which affects vibron transitions along 
the same spine. As a result, above a critical temperature 
I Jeff{i) | becomes the dominant hopping constant. In 
others words, the dressing effect strongly reduces the vi- 
brational exchange between different spines so that the 
small polaron tends to propagates almost along a sin- 
gle spine of hydrogen-bonded peptide units. Therefore, 
at biological temperature, the small polaron dynamics is 
rather similar to the dynamics described by the standard 
ID Davydov model. Note that the crossover between 
the two regimes, i.e. inter-spine and intra-spines vibron 
hops, is also induced by both an increase of the small po- 
laron binding energy and a modification of the #o angle 
which defines the helix backbone geometry. In this latter 
case, when 6>o tends to 2w/3, the vibrational exchanges 
between different spines vanish and only | J e //(3) | re- 
mains finite (see Fig. 5c). 

These features essentially originate from the coupling 
between the vibron and the acoustical phonons of the 
helix. Indeed, our calculations have revealed that the 
contribution to the dressing of the optical phonons is 
rather small, i.e. typically about two orders of magnitude 
smaller than the contribution of the acoustical phonons. 
This effect comes from the fact that the coupling constant 
S (n,T) is proportional to | Xi/Ki \ 2 whereas S a {n,T) 
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varies according to | X3/K3 \ 2 ■ Since xi 1S smaller than 
X3 and Ki is greater than K 3 , it is straightforward to 
show that S (n,T) « S a (n,T). 

At this step, the role played by the 6 angle as well 
as the fact that inter-spine or intra-spine vibron hops 
are subjected to a different dressing suggest that the ob- 
served features originate in the singular behavior of the 
acoustical phonons near the wave vector k = 2ir/3. In- 
deed, to understand more clearly this effect, let us assume 
an helix conformation with 80 = 2tt/3. As mentioned 
previously, in that case, the local frames n, n ± 3, n ± 6, 
... are equivalent. When k = 2n/3, the residues n, n±3, 
n±6, ... move in phase according to a zero frequency uni- 
form translation. As a consequence, the frequency of the 
nth amide-I vibration does not experience any modula- 
tion due to the motion of these residues. By contrast, the 
coupling constant S a (n, t) with residues located in differ- 
ent spines diverges since the phonon frequency is equal to 
zero. It results from this behavior around k = 2n/3 that 
the hopping constant along a given spine is slightly af- 
fected by the dressing effect whereas hopping constants 
between different spines vanish. When 9q decreases to 
100°, the equivalence between the local frames 11, n ± 3, 
n ± 6, ... is broken. Nevertheless, the motion of the 
residues along a given spine are more in phase than the 
motion due to residues located onto different spines. As 
a consequence, the modulation of a particular amide-I 
vibration due to the motion of residues located on differ- 
ent spines is more important than the modulation pro- 
duced by residues located in the same spine. Therefore, 
the dressing effect experienced by a vibron moving along 
a given spine is always smaller than the dressing effect 
which affects vibron hops between different spines. Note 
that the previous results corroborate the calculations per- 
formed by Hennig |l7j who have shown, within the soliton 
approach, that the preferred transport path takes place 
along a given spine whereas inter-spine energy transfer is 
suppressed. 

Finally, let us mention that the creation of a vibron 
in a 3D helix is accompanied by a deformation of the 
helix backbone different than the distortion which oc- 
curs in the ID Davydov model. Indeed, we have shown 
that the residues located on the spine where the vibron 
has been excited experience a more intense deformation 
than the other residues. This feature originates in the 
fact that inter-spine covalent bonds are more rigid than 
intra-spine hydrogen bonds. As a consequence, a distor- 
tion propagates more easily along a spine of hydrogen- 
bonded peptide units than along a covalent bond. The 
main consequence is that residues located on the excited 
spine do not react as residues located on the two other 
spines so that an asymmetric deformation of the helix 
backbone occurs. Due to this asymmetric behavior, the 
radial displacement of the residues located on the excited 
spine exhibits an important decrease whereas the radial 
displacement of the other residues just slightly increases. 
In the same way, both the longitudinal and the torsional 
deformations appear as kink-like solitons characterized 



by oscillations which account for the fact that the distor- 
tion essentially affects the excited spine, only. 

In that context, we basically recover the results ob- 
tained by Hennig 01 although our 3D model slightly 
differs from its cylinder model. Indeed, as Hennig, we 
have shown that the longitudinal distortion in a 3D helix 
is smaller than the distortion occurring in the standard 
ID model. This feature is due to the fact that in a 3D 
helix inter-spine covalent bonds tend to prevent the defor- 
mation of the hydrogen bonds. In addition, Hennig has 
shown that the helix radius decreases in the excited re- 
gion according to a deformation two orders of magnitude 
smaller than the longitudinal distortion. As explained 
previously, the helix model we use does not reduce to a 
cylinder. Therefore, when a vibron is excited on given 
site, the three spines react in an asymmetric way. Never- 
theless, the decrease of the helix radius was obtained, but 
for the residues belonging to the excited spine, only. In 
addition, we have shown that the amplitude of the radial 
deformation is about the same order of magnitude than 
the longitudinal distortion. In a marked contrast with 
Hennig who has assumed that the decrease of the radius 
originates in the optical phonons only, our model sug- 
gests that both optical and acoustical phonons strongly 
affect the radial deformation. In fact, our calculations 
have clearly shown that both the radial and the longitudi- 
nal distortions are essentially governed by the acoustical 
phonons so that they have the same order of magnitude. 



VI. CONCLUSION 

In the present paper, the vibron dynamics associated 
to amide-I modes in a 3D a-helix has been described. The 
helix was modeled by three spines of hydrogen-bonded 
peptide units linked via covalent bonds. To remove the 
strong coupling between the vibrons and the phonons as- 
sociated to the external motion of the residues, a Lang- 
Firsov transformation was applied and a mean field pro- 
cedure was performed to obtain the dressed vibron point 
of view. It has been shown that the vibron dynamics, 
due to the 3D nature of the helix, results from the com- 
petition between two kinds of hopping processes. The 
first kind refers to vibrational transitions between differ- 
ent spines whereas the second kind involves vibron hops 
along the same spine. Our study has revealed that several 
parameters such as the temperature, the small polaron 
binding energy and the helix backbone conformation, al- 
low for a transition between two regimes. Indeed, at 
low temperature or weak small polaron binding energy, 
the polaron behaves as an undressed vibron delocalized 
between the different spines. By contrast, at biological 
temperature or strong small polaron binding energy, the 
dressing effect strongly reduces the vibrational exchanges 
between different spines so that the polaron tends to 
propagate along a single spine, only. The occurrence of 
these two regimes also depends on the helix backbone 
geometry via the value of the 8q angle and it has been 
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shown that when 9q tends to 27r/3 inter-spine vibron hops 
are suppressed. Although the phonon spectrum exhibits 
both an acoustic branch and an optic branch, it has been 
shown that the previous features originate in the cou- 
pling between the vibrons and the acoustical phonons 
which exhibit a singular behavior near k = 2n/3. 

Finally, the previous results clearly show that the stan- 
dard ID Davydov model is a rather good approximation 
at biological temperature for which the dressing effect 
drastically reduces inter-spine vibron hops. In that case, 
for the vibron dynamics, the helix can thus be viewed as 
formed by three independent spines. 

To conclude, let us mention that the second paper of 
this series ^[ is devoted to the generalization of the 
present work when two vibrons are excited. In that case, 
both the intramolecular anharmonicity and the strong 
vibron-phonon coupling act as nonlinear sources which 
favor the occurrence of specific states called two- vibron 
bound states. As shown in paper II, these states exhibit 
an experimental signature within nonlinear pump-probe 
spectroscopy and they are expected to play a key role for 
energy storage in helices due to their formal resemblance 
with quantum breathers. 



where V'(nn') and V"(nn') denote the first and the sec- 
ond derivative of the pair potential between the residues 
n and n' , respectively. Note that the force constant ten- 
sor satisfies the well-known relation , § a p(nn') = 0. 



Within the representation of the displacements of the 
residues in the set of local frames, the transformed force 
constant tensor \I/ is expressed in terms of the local ro- 
tational matrices T(n) (Eq.JSJl), as 
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APPENDIX A: GENERAL EXPRESSION OF THE 
DYNAMICAL MATRIX 



In a general way, the force constant tensor Eq.10} is 
expressed as 
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where the matrix A(nn'), which depends on the index 
difference n — nl , is written in terms of the reduced pa- 
rameter h — h/Ro as 
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At this step, from Eq.(A3) and by using the notations 
introduced by Christiansen et al. |2Cj . the dynamical 
matrix D(k) = ^ \&(0n) exp(ikn)/M is finally defined 

as 



where 



Cll 



Cl2 



2a "[ 1 - cos(n6> )] 2 [l + cos(nk)] 

n 

2/3 n [l — cos(n^o) cos(nfc)] 

2a n [1 — cos n#o] sin(ro#o) sin(nfc) 



Cl3 



2/3 n sin(n^o) sin(nfc) 

2a n nh[l — cos uOq] sin(nk) 



(A5) 



D(k) 



Cll IC12 «Ci 3 ' 
-iCi2 C 2 2 C 2 3 
-iCl3 C 2 3 C33 



(A4) 



C22 = 2a n sin 2 (n9 ) [1 — cos nk] 



C23 



2/3 n [l — cos(n#o) cos(nfc)] 

2a n nh sin(n(9o)[l — cos nk] 



C33 = 2a n {nh) 2 [\ — cosnfc] 

n 

_l_ OR M _ nnalmhW 
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where the summation is performed over n 
and where 



1,2,3,. 



R? t 



MR(0n) 2 
V (On) 
MR(0n) 



V" (On) 



V(0n) 



R(0n) 



(A6) 



Note that, to reduce the number of parameters, we have 
assumed in our work that the helix conformation is sta- 
bilized when all the pair potentials are at equilibrium so 
that V'(On) = and K\ n \ = V"{0n). 



APPENDIX B: DISTORTION OF THE HELIX 
BACKBONE 

By disregarding the remaining coupling Hamiltonian 
AH, the lattice distortion onto the nth residue is ex- 
pressed, in the local frame, as 



= (^lUv^U- 1 ^) 



(Bl) 



where = 6^ |0) ® \ph) denotes the state associated 
to the creation of a polaron on the noth residue from the 
vacuum state \0) and where \ph) stands for the phonon 
wave function. By using Eq. i|13[) and after some straight- 
forward calculations this distortion is finally written as 



E 

kcr 



2NM cf/ 2 

kcr 



where cc. denotes the complex conjugate. From 
Eq.(B2), the z deformation corresponds to 6z n = x zn 
whereas the radius deformation 8R n and the angle defor- 
mation S9 n are defined as 



SR n — \ (Ro 



i) 2 +x 



On 



Ro 



(B3) 



_ x (Ro +x rn )smn0 o + x en cosn6 

56»„ = tan '(t^— : 2 = T> ~ n9 ° 

(Ro + x rn ) cos nt>o — xg n sm ntlo 
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